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ABSTRACT 

O 

O ■ We present sensitive phase-referenced VLBI results on the radio continuum 

emission from the z = 4.7 double source BRI 1202-0725. The observations were 



carried out at 1425 MHz using the Very Long Baseline Array (VLBA), the phased 
Very Large Array (VLA), and the Green Bank Telescope (GBT). 
^ ■ Our sensitive VLBI images of BRI 1202-0725 at 0.25 x 0.14 arcsec resolution 

show a continuum structure in each of its two components. Fitting Gaussian 



models to these continuum structures yield total flux densities of 315 ± 38 and 
250 ± 39 /uJy, for the northern and the southern components, respectively. The 
estimated intrinsic brightness temperatures of these continuum structures are 
~ 2 X 10^ K. Neither component is detected at the full VLBI resolution (29 mas 
x7 mas), with a 4cr point source upper limit of 40 fiJj beam~^, or an upper 
limit to the intrinsic brightness temperature of 6.7 x 10^ K. The highest angular 
resolution with at least a 4cr detection is ~85 mas. At this resolution, the images 
reveal a single continuum feature in the northern component of BRI 1202-0725, 
and two continuum features in the southern component, separated by 320 mas. 
This is similar to the structures seen in the high resolution images of the CO 
emission. 

The extent of the observed continuum sources at 1.4 GHz and the derived 
brightness temperatures are consistent with nuclear starbursts. Moreover, the 
absence of any compact high-brightness temperature source suggests that there 
is no radio-loud AGN in BRI 1202-0725. 
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1. INTRODUCTION 

Optical surveys such as the Sloan Digital Sky Survey (SDSS; York et al. 2000) and 
the Digitized Palomar Sky Survey (Djorgovski et al. 1999) have revealed large samples of 
quasi-stellar objects out to 2; ~ 6. Studies by Fan et al. (2002) have shown that at such 
a high redshift we are approaching the epoch of reionization, the edge of the "dark ages" , 
when the first stars and massive black holes were formed. 

Observations of high redshift QSOs at mm and sub-mm wavelengths have shown that 
a significant fraction (~ 30%) of the sources are strong emitters of far-IR radiation, with 

luminosities LpiR > 10"^^ Lq, and molecular gas masses greater than 10^° Mq (Brown & 
Vanden Bout 1991; Barvainis et al. 1994; Ohta et al. 1996; Omont ct al. 1996b; Guillotcau et 
al. 1997; Fraycr et al. 1999; Omont et al. 2001; Carilh et al. 2002, 2004). Moreover, the large 
reservoirs of warm gas and dust in these objects have led to the hypothesis that these are 
starburst galaxies with massive star formation rates on the order of 1000 Mq yr~^ (Hughes 
& Dunlop 1999; Lilly et al. 1999). 

An important question regarding these high-^ QSOs is whether the dominant dust 
heating mechanism is an AGN or a starburst. The high resolution of Very Long Baseline 
Interferometry (VLBI) observations permits a detailed look at the physical structures in the 
most distant cosmic sources. Also, sensitive VLBI continuum observations can be used to 
determine the nature of the energy source(s) in these galaxies at radio frequencies. 

To date, several high redshift QSOs have been imaged at milliarcsecond resolution (Prey 
et al. 1997, 2003; Momjian et al. 2003; Beelen et al. 2004). In this paper, we present 
sensitive VLBI observations of the z = 4.7 QSO BRI 1202-0725. This object was identified 
in the Automatic Plate Measuring survey (Irwin, McMahon, & Hazard 1991), and its optical 
spectra show strong associated Lya absorption (Storrie-Lombardi et al. 1996). BRI 1202- 
0725 was the first high redshift QSO in which thermal dust emission was detected (McMahon 
et al. 1994), 

Continuum emission at centimeter and millimeter wavelengths, as well as CO spectral 
line observations, showed that BRI 1202-0725 is a double source with an angular separation 
of 4". The total far-IR luminosity of this galaxy is LpiR = 4.2 x 10^^ Lq, with a dust 
temperature of ~ 70 K (Ohta et al. 1996; Omont et al. 1996a,b; Leech, Metcalfe, & Altieri 
2001; Carilh et al. 2002). 
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Various studies suggest that the observed double source may indicate a pair of interact- 
ing objects separated by 28 kpc, and that the large apparent luminosity of BRI 1202-0725 
may not be due to gravitational lensing (Omont et al. 1996b; Yun et al. 2000; Carilli et al. 
2002). Supporting this argument is the clear difference in the CO line widths and redshifts of 
the northern and southern components (Omont et al. 1996b; Carilli et al. 2002). Moreover, 
optical imaging of BRI 1202-0725 show a single point-like QSO with Mb = —28.5, with 
faint near-IR and Lya emission extending about 2'.'4 north of the QSO (Hu, McMahon & 
Egami 1996). The optical QSO is associated with the southern component seen in CO line, 
centimeter, and miUimeter, observations (CariUi et al. 2002). Most recently, Klamer et al. 
(2004) have proposed that the 4" double structure in BRI 1202-0725 could be due to jet 
induced star formation, with the northern component corresponding to a radio 'hot spot' 
(i. e. jet-shock or terminus). 

The CO (2 - 1) emission results with the VLA at 0'.'3 x 0''21 resolution (CariUi et al. 
2002) reveal a marginally detected source in the northern component of BRI 1202-0725 with 
an extent of > 0''5. The southern component of BRI 1202-0725 appears as two unresolved 
sources separated by 0'.'3 in these CO (2 — 1) observations. 

Throughout this paper, we assume a flat cosmological model with fi^ = 0.3, Qa = 0.7, 
and Hq — 65 km Mpc~^. In this model, at the distance of BRI 1202-0725, 1 milliarc- 
second (mas) corresponds to 7 pc. We also define the spectral index a, as a function of 
frequency {u), to be 5*1, oc i/". 

2. OBSERVATIONS AND DATA REDUCTION 

The VLBI observations were carried out at 1.4 GHz on 2003 April 16 and 19, using the 
Very Long Baseline Array (VLBA), the phased Very Large Array (VLA), and the Green 
Bank Telescope (GBT) of the NRAO^. Four adjacent 8 MHz baseband channel pairs were 
used in the observations, both with right and left-hand circular polarizations, and sampled 
at two bits. The data were correlated at the VLBA correlator in Socorro, NM, with 2 s 
correlator integration time and two independent correlation passes for the two components 
of BRI 1202-0725, which are separated by 4". The total observing time was 14 hr. Table 1 
summarizes the parameters of these observations. 

The observations employed nodding-style phase referencing, using the calibrator J 1202- 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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0528, with a cycle time of 4 min, 3 min on the target source and 1 min on the cahbrator. A 
number of test cycles were also included to monitor the coherence of the phase referencing. 
These tests involved switching between two calibrators, the phase calibrator J1202-0528 and 
the phase-check calibrator J1152-0841, using a similar cycle time to that used for the target 
source. 

The accuracy of the phase calibrator position is important in phase-referencing obser- 
vations (Walker 1999), as this determines the accuracy of the absolute position of the target 
source and any associated components. Phase referencing, as used here, is known to preserve 
absolute astrometric positions to better than ±0.01" (Fomalont 1999). 

Data reduction and analysis were performed using the Astronomical Image Processing 
System (AlPS) and Astronomical Information Processing System (AIPS++) of the NRAO. 
After applying a priori flagging, amplitude calibration was performed using measurements 
of the antenna gain and system temperature for each station. Ionospheric corrections were 
applied using the AIPS task "TECOR". The phase calibrator J1202-0528 was self-calibrated 
in both phase and amplitude and imaged in an iterative cycle. 

Images of the phase-check calibrator, J1152-0841, were deconvolved using two different 
approaches: (a) by applying the phase and the amplitude self-calibration solutions of the 
phase reference source J1202-0528 (Figure la), and (b) by self calibrating Jl 152-0841 itself, 
in both phase and amphtude (Figure lb). The peak surface brightness ratio of the final 
images from the two approaches gives a measure of the effect of residual phase errors after 
phase referencing (i.e. 'the coherence' due to phase referencing). At all times, the coherence 
was found to be better than 97%. 

The self-cahbration solutions of the phase calibrator, J1202-0528, were applied to the 
two correlation passes of the target source, BRI 1202-0725. Each of the two components of 
the target source was then deconvolved and imaged at various spatial resolutions by tapering 
the visibility data. 

In addition to these VLBI observations, 1.4 GHz VLA A-array (resolution ~ 1.5") 
observations carried out in August 1999 and November 2004 were reduced and analyzed for 
comparison purposes. 

3. RESULTS & ANALYSIS 

Figures 2 and 3 are naturally weighted images of the northern and southern components 
of BRI 1202-0725 at 1.4 GHz with a resolution of 252 x 141 mas (P. A.= -53°). The images 
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were obtained by applying a two-dimensional Gaussian taper falling to 30% at 1 MA in both 
the u- and ^^directions of the visibility data. The rms noises in these images are 28.3 and 
26.8 //Jy beam~^ for the northern and southern components, respectively. 

Table 2 shows the results of fitting Gaussian models to the observed source spatial pro- 
files using the AIPS-F-F task "IMAGEFITTER" . The source names, redshifts, and positions 

are listed in columns 1, 2, 3, and 4, respectively. Columns 5 and 6 show the peak and total 
flux densities of the fitted Gaussian functions. Columns 7, 8, and 9 list the fitted half-power 
ellipse axes and the position angles. The intrinsic brightness temperatures (corresponding 
to a rest-frame frequency of 8 GHz) of the radio continuum structures seen in the northern 
and southern components of BRI 1202-0725 are (2.2 ± 0.3) x 10^ and (1.7 ± 0.3) x 10^ K, 
respectively. The total radio flux density of BRI 1202-0725 at 1.4 GHz, i. e., the sum of 
the flux densities of both components, is 565 ± 55 /xJy, as measured with our sensitive VLBI 
array. 

No continuum emission is detected at the the full resolution of the array, which is 
29 X 7 mas (PA=— 6°), indicating the absence of any radio continuum emission with flux 
densities of > 4(7 ~ 40 //Jy beam~^. This implies an upper limit to the intrinsic brightness 
temperature of 6.7 x 10^ K for any compact radio source in BRI 1202-0725. We emphasize 

again that our coherence tests during these observations using two VLBI calibrators show 
that the lack of a strong point source in BRI 1202-0725 at the full resolution of the array 
cannot be due to the phase referencing procedure. 

Figures 4 and 5 show the two components in BRI 1202-0725 at the highest angu- 
lar resolutions for which there are at least 4cr detections. These detections correspond to 
resolutions of 81 x 57 mas in position angle —11° for the northern component (Figure 4; 
a = 25.7 /jJj beam^^), and 111 x 83 mas in position angle —37° for the southern component 
(Figure 5; a = 25.5 ^Jy beam~^). These images were obtained by applying two-dimensional 
Gaussian tapers falling to 30% at 2.6 and 1.8 MA in both the u- and ^^directions of the 
visibility data, for the northern and southern components, respectively. 

The intermediate angular resolution image of the northern component (Figure 4) shows a 
single continuum structure at 4.3o" level. However, the southern component (Figure 5) seems 
to be resolved into two continuum structures separated by 320 mas. The stronger continuum 
source in Figure 5 is detected at 4.8(7 level, and the weaker at 3.7a. The derived intrinsic 
brightness temperature of the continuum structure in Figure 4 (northern component) is 
(7.5 ±1.9) X 10^, and of the structures in Figure 5 (southern component) are (4.5 ±1.0) x 10^ 
and (3.5 ± 1.0) x 10^ K. The extent of these cont imnim sources range between 0.4 and 1 kpc 
on the plane of the sky. The single, resolved structure of the northern component, and the 
double structure of the southern component seen in these VLBI images are similar to the 
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structures seen in the CO (2 — 1) emission observations of BRI 1202-0725 (Carilli et al. 
2002). In particular, the CO (2 — 1) emission from the southern component shows a 300 
mas separation double morphology along the same position angle as that seen in the radio 
continuum image of Figure 5. 

In addition to the VLBI results, we analyzed VLA A-array data from two epochs, 
August 1999 and November 2004. The VLA flux densities, which are listed in Table 3, are 
comparable to the VLBI measurements reported in Table 2. 

4. DISCUSSION 

At a moderate angular resolution (252 x 141 mas; Figures 2 & 3), the two components 
of the high z QSO BRI 1202-0725 show radio continuum emission at the observed 1.4 GHz, 
which corresponds to a rest frequency of 8 GHz at 2; = 4.7. The northern and southern 
components constitute 56 and 44% of the total flux density of BRI 1202-0725, respectively. 
The intrinsic brightness temperature values of these continuum structures are on the order 
of 2 X 10^ K. 

At the full resolution of our array (29 x 7 mas), the two components of BRI 1202- 
0725 are resolved out and do not show any single dominant source of very high brightness 
temperature (< 6.7 x 10^ K). This is in contrast to the results obtained by Momjian, Petric, 
& Carilli (2004) on a sample of three high- 2; radio-loud quasars which were imaged with the 
VLBA. In each of these 2; > 4 quasars, a radio-loud AGN dominates the emission at 1.4 GHz 
on a few mas size scale, with intrinsic brightness temperatures in excess of 10^ K. 

Condon et al. (1991) have derived an empirical upper limit to the brightness temperature 
for nuclear starbursts of order 10^ K at 8 GHz, while typical radio-loud AGN have brightness 
temperatures exceeding this value by at least two orders of magnitude. These authors also 
present a possible physical model for this limit involving a mixed non-thermal and thermal 
radio emitting (and absorbing) plasma, constrained by the radio-to-FIR correlation for star- 
forming galaxies. 

For the two components of BRI 1202-0725, the derived intrinsic brightness temperatures 
from our VLBI observations are typical of starburst galaxies. Also, the hnear extents of the 
radio continuum emission regions in BRI 1202-0725 are typical of local starburst powered 

Ultra-Luminous IR galaxies (ULIRGs; Sanders & Mirabel 1996), such as Arp 220, Mrk 273, 
and IRAS 17208-0014 (Smith et al. 1998; Carilh & Taylor 2000; Momjian et al. 2003). 
Moreover, the radio continuum morphology of each component in BRI 1202-0725 is similar 
to that seen in high resolution images of the CO emission (Carilli et al. 2002). Co-spatiality 
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of the radio continuum emission with the molecular gas (and presumably dust) is expected 
for star forming galaxies (CariUi et al. 2003). 

The physical characteristics of the QSO BRI 1202-0725 are also similar to one of the 
most luminous IR sources, J1409+5628, located Bit z — 2.58 (Beelen et al. 2004). High 
angular resolution observations at 1.4 GHz showed that the radio continuum emission in 
this optically very bright QSO is dominated by an extreme nuclear star burst, with a massive 
star formation rate of (3 — 8) x 10^ Mq yr~^. 

Another argument in favor of a massive starburst origin for the radio continuum emission 
from BRI 1202-0725 is the spectral index between the 1.4 and 350 GHz {al^2) for the overall 

flux density of its two components. This spectral index, af^^, corresponds to that measured 
between the non-thermal synchrotron radiation at 1.4 GHz and the thermal emission from 
warm dust at 350 GHz. For BRI 1202-0725, afj" = +0.96 ±0.04, and the predicted value for 
a starburst galaxy at z = 4.7 is aff = +1-03 ± 0.10 (Yun et al. 2000; Carilh & Yun 2000). 

Yun et al. (2000) suggested possible flux density variation at 1.4 GHz for BRI 1202- 
0725. Such variation would be hard to reconcile physically with the kpc-scale spatial extent 
observed with VLBI for the radio components. However, the variability conclusion was based 
on comparison of VLA observations made with radically different resolutions (A through D 
array; or 1" and 60"). We have reexamined the variabihty using observations made with the 
VLA A-array in August 1999 and November 2004. The total flux densities of each component 
in BRI 1202-0725, as reported in Section 3, do not show any sign of variabihty during a five 
year period. Moreover, the measured fiux densities with the VLA A-array are consistent 
with our VLBI results (Table 2), suggesting that the radio continuum emission at 1.4 GHz 
is confined within the extent of the structures seen in the VLBI images. 

While the measured 1.4 GHz radio fiux densities in BRI 1202-0725 are consistent with 
extreme nuclear starburst activity, detecting individual supernovae (SNe) at z = 4.7 is 
unlikely. Smith et al. (1998) reported the detection of several luminous radio SNe in the 
prototype ULIRG Arp 220 with flux densities between 0.2 and 1.17 mJy and angular sizes on 
the order of 0.25 mas. At the distance of BRI 1202-0725, the flux densities of such luminous 
radio SNe would be between (0.7 — 4) x 10~^ /xJy. These values are four orders of magnitude 
lower that the rms noise levels achieved in our VLBI observations. 

In general, the single source morphology of the northern component and the double 
source morphology of the southern component in BRI 1202-0725, and their starburst nature, 
lead us to conclude that the northern component could be similar in nature to single nucleus 
ULIRGs, such as IRAS 17208-0014 (Momjian et al. 2003), and the southern component to 
double nucleus ULIRGs, such as Arp 220 (Smith et al. 1998). Moreover, these two extreme 
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starburst galaxies (i. e. the two components of BRI 1202-0725), that are separated by 4", 
may actually be interacting as suggested by various studies (Omont et al. 1996b; Yun et al. 
2000; CariUi et al. 2002). 
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Fig. 1. — Continuum images of the phase-check cahbrator J1152-0841 at 1.4 GHz: a) ob- 
tained by applying the phase and the amphtude self-cahbration solutions of the phase ref- 
erence source J1202-0528, b) obtained by self cahbrating Jl 152-0841 itself, in both phase 
and amplitude. The restoring beam size in both images is 17.2 x 5.8 mas in position angle 
— 1°. The contour levels are at —3, 3, 6, 12, . . ., 192 times the rms noise level in the phase- 
rcfcrcnccd image (top), which is 1.2 mJy beam~^. The gray-scale range is indicated by the 
step wedge at the right side of each image. 
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Fig. 2. — Naturally weighted 1.4 GHz continuum image of the northern component of 
BRI 1202-0725 at 252 x 141 mas resolution (P. A.=-53°). The peak flux density is 
232 iiJy beam~^, and the contour levels are at —3, 3, 4, 5, 6, 7, 8 times the rms noise 
level, which is 28.3 //Jy beam~^. The gray-scale range is indicated by the step wedge at 
the right side of the image. The reference point (0, 0) is «(J2000.0) = 12^05™23?9785, 
(5(J2000.0) = -07°42'29!'7758. A two dimensional Gaussian taper falhng to 30% at 1 MA 
was applied. 
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Fig. 3. — Naturally weighted 1.4 GHz continuum image of the southern component of 
BRI 1202-0725 at 252 x 141 mas resolution (P. A.=-53°). The peak flux density is 
177 /iJy beam~^, and the contour levels are at —3, 3, 4, 5, 6 times the rms noise level, 
which is 26.8 /iJy beam"^. The gray-scale range is indicated by the step wedge at the 
right side of the image. The reference point (0, 0) is a(J2000.0) = 12^05°'23?1187, 
(5(J2000.0) = -07°42'33!'l391. A two dimensional Gaussian taper falling to 30% at 1 MA 
was applied. 
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Fig. 4. — Naturally weighted 1.4 GHz continuum image of the northern component 
of BRI 1202-0725 at 81 x 57 mas resolution (P. A.=-ll°). The peak flux density is 
108 fxJy beam~^, and the contour levels are at —3, 3, 4 times the rms noise level, which 
is 25.7 yuJy beam~^. The gray-scale range is indicated by the step wedge at the right side 
of the image. The reference point (0, 0) is a(J2000.0) = 12^05"^23?9785, 5(J2000.0) = 
— 07°42'29'.'7758. A two dimensional Gaussian taper faUing to 30% at 2.6 MA was applied. 
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Fig. 5. — Naturally weighted 1.4 GHz continuum image of the southern component of 
BRI 1202-0725 at 111 x 83 mas resolution (P. A.=-37°). The peak flux density is 
110 fxJy beam~^, and the contour levels are at —3, 3, 4 times the rms noise level, which 
is 25.5 /xJy beam~^. The gray-scale range is indicated by the step wedge at the right side 
of the image. The reference point (0, 0) is a(J2000.0) = 12^05"^23!1187, 5(J2000.0) = 
— 07°42'33'.'l391. A two dimensional Gaussian taper faUing to 30% at 1.8 MA was applied. 
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Table 1. PARAMETERS OF THE VLBI OBSERVATIONS OF BRI 1202-0725 



Parameters Values 



Observing Date 2003 Apr 16 & 19 

Total observing time (lir) 14 

Phase ealibrator J1202-0528 

Phase-refereneing cycle time (min) 4 

Frequency (MHz) 1425 

Total bandwidth (MHz) 32 



Table 2. SOURCE PARAMETERS in FIGURES 2 & 3 



Gaussian Component Parameters 



Source z R.A. (J2000) Decl. (J2000) Peak Total Major Axis''' Minor Axis'^' P. A.<^ 

(h m s) (° ' ") (/iJy beam-l) (//Jy) (mas) (mas) (°) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



Northern 4.6916 12 05 22.977 -07 42 29.750 234 ± 28 315 ± 38 288.6 166.5 115 

Southern 4.6947 12 05 23.1167 -07 42 33.109 176 ± 27 250 ± 39 297.4 170.6 130 

''At half maximum. 

''Gaussian fit error is < 0.3 mas. 

^Gaussian fit error is < 0.001°. 
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Table 3. MEASURED 1.4 GHz FLUX DENSITIES WITH THE VIA A- ARRAY 



Source 


Date 


Si. 4 GHz 






(A'Jy) 


Northern 


Aug. 1999 


365 ± 27 




Nov. 2004 


282 ± 48 


Southern 


Aug. 1999 


230 ± 34 




Nov. 2004 


203 ± 26 



